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Abstract 

The a-decay chain of ^^^118, first proposed in the Berkeley cold fusion 
experiment ^°'^Pb(^^Kr,ln) and now retracted, is calculated by using the pre- 
formed cluster model (PCM) of one of us (RKG). Also, the possible branch- 
ings of a-particles to heavier cluster decays of all the parents in this chain 
are calculated for the first time. The calculated Q-values, penetrabilities 
and preformation factors for a-decays suggest that the ^^^114 nucleus with 
Z=114, N=171 is a magic nucleus, either due to the magicity of Z=114, or of 
N=172 or of both. The N=172 is proposed to be a magic number in certain 
relativistic mean-field calculations, but with Z=120. The calculated cluster 
decays point to new interesting possibilities of ^^C decay of the ^^^112 par- 
ent, giving rise to a (reasonably) deformed Z=106, N=161, ^^''106 daughter 
(N=162 being now established as the deformed magic shell) or to a doubly 
magic ''^Ca cluster emitted from any of the parent nucleus in the a-decay 
chain. Apparently, these are exciting new directions for future experiments. 
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The synthesis of Z=118 element in the cold fusion reaction ^'^^Pb(^^Kr,ln) 
via the observed a-decay chains had created much excitement recently. This 
reaction was first made at Berkeley [1], establishing three decay chains, with 
a resulting very high fusion cross section of 2.2lo;8 pb compared to the lim- 
iting value of ~ 1 pb for the cold fusion reactions leading to the heaviest 
Z=112 element. However, this experiment is now retracted [2, 3], since 
many other subsequent attempts [4, 5, 6, 7] at various other laboratories 
(GSI,RIKEN,GANIL) around the World failed to reproduce these data and 
"one event" upper limit of <0.5 pb has now been put for this reaction. Such 
a resulting situation has some consequences for what has been the cause for 
a large excitement for nuclear structure studies in the region of superheavy 
elements [8, 9, 10, 11, 12, 13, 14], discussed below. 

The measured large cross section for Z=118 element, now retracted, was 
considered [13] as a possible signature of our approaching the centre of island 
of stabihty for superheavy elements (SHE) around Z=120, predicted by the 
well founded relativistic mean-field (RMF) calculations [8, 9, 10, 11, 12]. 
However, the lowering down of the fusion cross section for this reaction to 
<0.5 pb means that we must go up the ladder of SHE rather steadily, as 
was proposed by another calculation of some of us [15], based on the well 
accepted Quantum Mechanical Fragmentation Theory [16, 17, 18, 19, 20, 21]. 
The surprises, if any, were expected to lie in the overshooting of this centre of 
island of SHE by means of (neutron-rich) radioactive nuclear beams for the 
magic N=184. Then the question would arise for protons, whether Z=110 
or 114 is magic, as has been predicted since the early days of this subject, 
or it is around Z=120, as is predicted more recently by the above mentioned 
RMF calculations. 

In this paper, we attempt to look for an answer to the question raised 
in the last paragraph above, regarding the nuclear structure effects in SHE, 
as well as to the cause for the failure of Z=118 experiment. We do this 
by analysing theoretically the a-decay chain for ^^^118. This is only an 
exploratory study and can be extended to other heavy nuclei in this region. 
Also, we have calculated for the first time the possible branching of a-decay 
to any heavy cluster decay, at any stage of the a-decay chain of ^^^118. 
Such a process of, so-called, cluster radioactivity should open new vistas for 
the decay studies of SHE, with a possible landing at some new or known 
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magic daughter nucleus. We have used here for our decay calculations the 
Preformed Cluster Model (PCM) of Gupta [22, 23, 24]. 

The preformed cluster model (PCM) uses the dynamical collective coor- 
dinates of mass and charge asymmetries 

77= (^i-A)/(A + A) 

and 

first introduced in the Quantum Mechanical Fragmentation Theory [16, 17, 
18]. These are in addition to the usual coordinates of relative separation 
R and deformations /3j {i = 1,2). Then, in the standard approximation of 
decoupled R- and ?7-motions [22, 23, 24, 25], the decay half-life 71/2 or the 
decay constant A in PCM is defined as 

A = 1^ = Po-oP. (1) 

Here Pq is the cluster (and daughter) preformation probabihty and P the 
barrier penetrability which refer, respectively, to the rj and R motions. The 
i^o is the barrier assault frequency. The Pq are the solutions of the stationary 
Schrodinger equation in rj, 

which on proper normalization are given as 

Po^^\^^'\v{A))n2/A), (3) 

with i=l or 2 and i/=0, 1,2,3.... We are interested here only in the ground 
state solution (I'—O) since the a (as well as the proposed heavy cluster) 
emissions in the considered decay chain are the ground state decays. Eq. 
(2) is solved at a fixed R — Ra — Ct{— Ci + C2) + d, the first turning 
point in the WKB integral for penetrability P (see Fig. 1 and Eq. 5), 
since this value of R (instead of i? = i?o, the compound nucleus radius) 
assimilates to a good extent the effects of both the deformations of two 
fragments and neck formation between them [26]. In this way, the two- 
centre nuclear shape is simulated through a neck-length parameter d added 
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to scission configuration, which for actinides is nearly zero [26] , and is taken 
to be so for superheavy nuclei. The role of deformation in the scattering 
potential V(R) is shown [26] to lower the interaction barriers but not the 
relative formation yields. The Cj are Siissmann central radii Cj = Ri—{1/Ri), 
with the radii Ri = 1.28^^/^ - 0.76 + 0.8^"^/^^- 

The fragmentation potential VR{ri) in (2) is calculated simply as the sum 
of the Coulomb interaction, the nuclear proximity potential [27] and the 
ground state binding energies of two nuclei, 

V{Ra, V, Vz) = - E B{A, Z,, A) + Vp, (4) 

i=i 

with B's taken from the 1995 experimental compilation of Audi and Wapstra 
[28] and from the 1995 calculations of MoUer et al. [29] whenever not available 
in [28]. Thus, full shell effects are contained in our calculations that come 
from the experimental and/or calculated [29] binding energies. The shell 
effects in the calculated binding energies of MoUer et al. [29] are obtained 
in Strutinsky way [30] by using the folded- Yukawa single-particle potential 
and macroscopic finite-range droplet model (FRDM). The model parameters 
are fitted to the ground-state masses of 1654 nuclei, ranging from ^^O to 
^^^106, and to 28 fission-barrier heights. Hence, their extrapolation to heavier 
elements, studied here in this paper, is expected to be a realistic one. Note 
that the familiar magic numbers N=Z=50, 82 and N=126 are given by these 
calculations and there is no assumption made about the magic numbers in 
the extrapolated region of SHE. The center of the superheavy region in these 
calculations is found to be located at ^'^'^llSiyg [29]. 

The charges Zi and Z2 in (4) are fixed by minimizing the potential in r]z 
coordinate, which automatically minimizes the /3i coordinates. Note that the 
minimized /3j's are not always for the spherical nuclei since the total binding 
energy Bi + B2 of the decay products is minimized and not their individual 
Bi or B2. The Coulomb and proximity potentials in (4) are for spherical 
nuclei. The mass parameters Bnri{r)), representing the kinetic energy part in 
(2), are the classical hydrodynamical masses [31]. The shell effects in masses 
are shown [24] not to affect the order of calculated preformation yields. 

The WKB tunnelling probability, calculated for the tunneling path shown 
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in Fig. 1, is P = PiPb with 



p. ^ exp[-- / {2i,[V{R) - V{Ri)]} 



2 fR^ 



(5) 



(6) 



J Ri 



These integrals are solved analytically [23] for R^, the second turning point, 
defined by V{Rb) = Q- value for the ground-state decay. 
The assault frequency fo in (1) is given simply as 



with E2 = {Ai/A)Q, the kinetic energy of the lighter fragment, for the Q- 
vahie shared between the two products as inverse of their masses. Eq. (7), 
used here, results in z/q ~ 3 x lO^^s^^, whereas the more often used value in 
literature is ~ lO^^s"^ for even parents and ~ lO'^^s^^ for odd parents [32]. 

Figure 2 shows the calculated (logarithms of) a-decay half-lives, logi^Ty^ 
(s), as a function of the masses of parent nuclei for the whole decay chain 
of ^^^118. Also shown in Fig. 2 are the results of another recent calculation 
by Royer [33] based on the generalized liquid drop model (GLDM) with 
the binding energies for Q-values taken from the Thomas-Fermi (TF) model 
[34]. The TF model uses the well known Seyler-Blanchard effective nucleon- 
nucleon interaction where the momentum-dependent and density-dependent 
terms are also included. The model parameters are fitted again to ground- 
state masses of 1654 nuclei with N,Z>8 and the 40 fission-barrier heights. 

We notice in Fig. 2 that the calculated numbers for both the models 
present an interesting result: the a-decay half-life for Z=114, A=285 is very 
high, which means that the parent nucleus ^^^114 is very stable against a- 
decay. This stability can be attributed to either the magicity of protons at 
Z=114 or of neutrons aX N 172 or to both, perhaps more so to Z—llA 
since N=172 is predicted to be magic only when Z=120 [8, 9, 10, 11, 12, 14]. 
The predicted half-lives in the two models differ by about four orders of 
magnitude, which is partly due to different choices of i^o-values. 

The stabihty at Z=114 in the calculations arises due to the Q-values 
involved, as is evident from Fig. 3(a) where the calculated Q-values are 
plotted as a function of the parent nuclear masses ^Z for both the PCM 
and GLDM models, as well as that of another calculation by Smolanczuk 
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[35] . The Q- values are similar for the PCM and GLDM models and are very 
small (minimum) for the a-decay of ^^^114 parent. However, this is contrary 
to the predictions of Smolahczuk [35], which guided the Berkeley retracted 
experiment. Smolahczuk [35] predicted the Q- values as an ever increasing 
function of the parent nucleus charge (or mass), and became an apparent 
cause for the failure of this experiment [1] . 

Smaller Q- value should also mean a relative decrease in the penetrability 
P. This is shown to be the case in our calculations presented in Fig. 3(b) 
where —IoqiqP vs. "^Z is plotted. It is further interesting to find that the 
calculated preformation factor Pq in Fig. 3(c) is also minimum for the a- 
nucleus preformation in ^^^114 parent. The Pq factors in present calculations 
are shown smaller by a few orders of magnitude (~ 10""^), compared to that 
for the actinides [36] . 

So far, the Z—llA element (A=287-289) is synthesized only in hot fusion 
reactions made at Dubna, using a ^^Ca beam on 242,244p^ targets [37, 38, 39] , 
which are found to result in larger production cross sections [4]. The cold 
synthesis of Z=114, in the proposed reaction '^^Ge +^°^ Pb [19], is still an 
open question experimentally where the measured cross section could throw 
some light on its magic structure, if any. Apparently, the cold identification 
of Z—llA will prove a testing ground for many structure calculations. 

Figures 4 and 5 give the results of our calculations for heavy cluster decays 
of each of the parents in the ct-decay chain of ^^^118, for a few illustrative 
clusters. The choice of clusters is based on the minima in the fragmentation 
potentials V{rj) and hence for the cases of largest preformation factors Pq, 
illustrated as an example for Z=118 in Fig. 6. The heaviest cluster included 
here in Figs. 4 and 5 is with Z—20 (^^~^°Ca) because the earher calculations 
on PCM show that the two processes of cold fission and cluster decay are 
almost indistinguishable for clusters heavier than of mass A2 > 48 [40] . 

First of all we notice in Fig. 4(a) that the Q- value increases as the size of 
cluster increases, but is almost independent of the parent mass (the increase 
with parent mass is smooth, linear and with a very small slope). However, 
the calculated cluster decay half-hves, ^^S- ^(b) present some 

interesting results: (i) The decay half-hves for all the hght clusters, other 
than the ^^C from Z=112-116 parents, are rather high and hence the studied 
parents could be said as stable against most of the hght cluster decays. The 
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shell stabihzing effect, if any, is seen for ^°Be decay of ^^^116 or ^"^108 nucleus, 
since show strong peaking at these two parent nuclei, (ii) The ^^C 
decay, in particular from ^^^112 parent, seems to present an interesting case 
of (possibly, a reasonably) deformed magic daughter ^^^106 with N ^ 162. 
This means to say that for the known deformed magic shell at N=162 [4], 
the ^^^Sg, not yet synthesized, could also be deformed and observed as ^^C 
decay of ^^^112 parent, (iii) The heavier clusters ^^~^°Ca are predicted to 
decay with further smaller half-lives and hence present themselves as further 
interesting cases of cluster decay measurements. The calculated half-lives 
for ^^"^'^Ca decays lie far below the present hmits of experiments, which go 
upto ~ lO^^s [32] for nuclei where enough atoms are available. Here the 
closed shell effects of cluster (not daughter) are playing the role, for which 
at present no measurements exist in radioactive cluster decay studies. The 
heaviest cluster observed so far is ^^Si from ^^^Pu parent. 

Figure 5 gives the cluster preformation (Pq) and penetration (P) proba- 
bihties. Knowing that T1/2 is a combined effect of both Pq and P {vq being 
almost constant), we notice in Fig. 5 that though ^'^Be is better preformed 
(larger Pq) than both ^^C and Ca nuclei, but its penetration probabilty P is 
so small that the 71/2 for ^°Be decay is much larger than for either of the 
two other clusters. Thus, in experiments one should consider the possibility 
of ^^C and/ or ^^Ca decays in addition to a-decay or fission of any of the 
parents in ^^^118 ct-decay chain. 

Summarizing, we have calculated the a-decay of ^^^118 and its subsequent 
parents ending the chain in ^^^106, as well as the heavy cluster decays of all 
the parents in the a-decay chain. Though the experimental data for this a- 
decay chain is retracted, the calculated a-decay half-lives are found to contain 
interesting nuclear structure information. For the Q- values calculated from 
experimental binding energies, supplemented by the Finite Range Droplet 
Model (FROM) or the Thomas-Fermi (TF) Model calculations, the a-decay 
half lives show that there is a magic shell structure at either or both Z=114, 
Nf=:il72, since the calculated T1/2 value for ct-decay of ^^^114 shows a strong 
peaking structure. This is supported by a small Q-value and small (deeper 
minima) penetrability and preformation factors for a in the ^^^114 parent. 

The cluster decay calculations for a decaying superheavy nucleus formed 
in heavy ion reactions is made for the first time, with a view to see if there 
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is any branching of a-decay to another hght nucleus due to the (spherical/ 
deformed) magicity of the corresponding heavy daughter nucleus or due to 
the magicity of the hght nucleus itself. Interesting enough, the ^^^118 decay 
chain offers two such possibilities: firstly, the ^^C decay of the inbetween par- 
ent ^^^112, and secondly, the doubly magic ^^Ca decay of any parent nucleus 
obtained after the Q;-decay(s) in the investigated chain. The first possibil- 
ity points out to the deformed magicity of the daughter product ^^^106 at 
N=162 and the second possibility to a first time observation of the doubly 
magic emitted cluster ^^Ca. These are exciting new possibilities for future 
studies. 
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Figure Captions 

Fig. 1: The scattering potential for a-decay of ^^^118, calculated as the 
sum of Coulomb and nuclear proximity potential. The tunneling path used 
by the PCM is also showm with the first and second turning point radii 
marked as Ra and Rb, respectively. 

Fig. 2: The a-decay half-hves calculated on PCM and compared with 
GLDM model [33] plotted as a function of the parent nucleus mass for the 
a-decay chain of ^^^118. 

Fig. 3: (a) The Q- values calculated on FRDM [29], compared with TF 
[34] and Smolahczuk [35] calculations, (b) the penetration probabilty P, and 
(c) the preformation factor Pq calculated on PCM, plotted as a function the 
parent nucleus mass for the ct-decay chain of ^^^118. 

Fig. 4: The calculated Q-values on FRDM [29] and decay half-hves on 
PCM for some cluster decays of the parents of a-decay chain for ^^^118 plot- 
ted as a function of the parent nucleus mass. 

Fig. 5: Same as for Fig. 4, but for the penetration probabilty P and the 
preformation factor Pq. 

Fig. 6: (a) The fragmentation potential V{rj) for ^^^118, using the exper- 
imental and FRDM binding energies [28, 29]. Only half of the potential as a 
function of the cluster mass A2 is shown. The other (symmetrical) half is a 
function of the daughter mass Ai. Some of the minima are marked, showing 
the cluster mass. 

(b) The preformation factor Pq, plotted as a function of both the cluster and 
daughter masses, for the decaying ^^^118 parent. 
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